! ******************************
/units,si 
/title, Fatigue analysis of cylinder with flat head
! ***************趨***************
Di=1000                                ! Ͳھ
t=20                                   ! Ͳ
hc=nint(4*sqrt(Di/2*t)/10)*10          ! ģͲ峤
tp=60                                  ! ƽͷ
r1=10                                  ! ƽͷԲ뾶
r2=10                                  ! ƽͷڲӦͷŲԲ뾶
exx=2e5                                ! ϵģ
mu=0.3                                 ! ϲɱ
p1=2                                   ! ߹ѹ
p3=2.88                                ! ˮѹѹ
n1=2e4                                 ! /ѹѭ
n2=5                                   ! ˮѹ
! ***************ǰ***************
/prep7
et,1,82                                ! 趨Ԫ
keyopt,1,3,1                           ! 趨ܶԳѡ
mp,ex,1,exx                            ! ϵģ
mp,nuxy,1,mu                           ! ϲɱ
!  *******  ģ  *******
k,1,0,0                                ! ؼ
k,2,Di/2+t,,
k,3,Di/2+t,-(tp+hc)
k,4,Di/2,-(tp+hc)
k,5,Di/2,-tp
k,6,Di/2-r2,-tp                        ! ӦͷŲԲĹؼ
k,7,0,-tp
l,1,2                                  ! 
l,2,3
l,3,4
l,4,5
l,5,7
l,7,1
LFILLT,1,2,r1                          ! Բ   
al,all                                 !  
CYL4, kx(6),ky(6), r2,180              ! ӦͷŲ   
ASBA,1,2                               ! 
wprot,,,90                             ! תƽ
wpoff,,,kx(6)-3*r2                     ! ƶƽ 
asbw,all                               ! ùƽи
wprot,,90                              ! תƽ
wpoff,,,tp+r2                          ! ƶƽ 
asbw,all                               ! ùƽи

esize,5                                ! 趨Ԫߴ
MSHKEY,1                               ! 趨ӳʷ
amesh,1                                ! ӳʷ1
amesh,3                                ! ӳʷ3
esize,2                                ! 趨Ԫߴ
MSHKEY,0                               ! 趨ʷ
amesh,4                                ! ʷ4
fini                                   ! ˳ǰ
! ******************************
/solu                                  ! Ͳ˲ʩԼ
dl,3,,uy                               ! Ͳ˲ʩԼ
dl,6,,symm                             ! ƽͷԳʩӶԳԼ

time,1                                 ! غɲ1
lsel,s,,,8                             ! ѡڱ߶
lsel,a,,,11,13
lsel,a,,,15
cm,lcom1,line                          ! ڱ
SFL,all,PRES,p1,                       ! ڱʩѹ
alls                                   ! ȫѡ
solve                                  ! 
fini                                   ! ˳
! ******************************
/post1                                 ! 
FTSIZE,1,2,2,                          ! 趨ƣλ¼غ

FP,1,1e1,2e1,5e1,1e2,2e2,5e2           ! ƣSN
FP,7,1e3,2e3,5e3,1e4,2e4,5e4
FP,13,1e5,2e5,5e5,1e6, ,
FP,19, ,
FP,21,4000,2828,1897,1414,1069,724  
FP,27,572,441,331,262,214,159   
FP,33,138,114,93.1,86.2, ,  
FP,39, ,

! ******  ˮѹѭ ******
fs,4760,1,1,1,0,0,0,0,0,0              ! ڵ4760ӦһغɵӦ
set,1,last                             ! һغɲ
FSNODE,4760,1,2                        ! ڵ4760ӦڶغɵӦ
fe,1,n2,p3/p1                          ! 趨¼ѭغɱϵ
! ******  /ѹѭ ******
fs,4760,2,1,1,0,0,0,0,0,0              ! ڵ4760ӦһغɵӦ
set,1,last                             ! һغɲ
FSNODE,4760,2,2                        ! ڵ4760ӦڶغɵӦ
FE,2,n1,1,                             ! 趨¼ѭغɱϵ 
FTCALC,1                               ! ƣͼ㣨¼ʹϵ
fini
